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With  the  current  emphasis  on  composite  materials,  it  has  become  necessary  to 
use  the  thick -wall  cylinder  equations  for  cylindrically  orthotropic  materials. 
This  report  is  a  preliminary  investigation  of  these  equations.  The  most 
important  result  is  that  many  of  the  old  ideas  of  high  pressure  cylinders  will 
have  to  be  changed.  When  the  cylinder  ie  very  thin  (wsll  ratio  near  1.0),  the 
simple  thin-wall  equations  are  adequate.  However,  aa  the  thickness  increases, 
stress  variation  through  the  wall  becomes  large  and  the  full  thick-wall 
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solution  is  necessary.  In  orthotropic  (composite)  cylinders  this  transition 
happens  at  a  lower  wall  ratio.  Furthermore,  the  maximum  useful  wall  ratio 
may  become  smaller  for  a  composite  cylinder.  -'K*  >  u  •  \ 
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tha  material  compliance  matrix 

the  outer  radius 

the  cylinder  geometry  parameter 

simplified  orthotropic  material  parameter 

engineering  modulus  in  the  hoop  direction 

engineering  modulus  in  the  radial  direction 

a  material  constant  for  .a!  load  equations 

matrix  row  index 

matrix  column  index 

orthotropic  material  parameter 

internal  pressure 

axial  force 

plane-strain  axial  force 
external  pressure 
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calculated  constant  for  the  axial  load  equations 
wall  ratio 
radial  direction 

circumferential  or  hoop  direction 
axial  direction 

compliance  for  the  cylindrical  problem 
strain  in  the  i  direction 


stress  in  the  i  direction 
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»  matrix  position  for  the  x  direction 
«  matrix  position  for  the  y  direction 
*  matrix  position  for  the  z  direction 
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[  ]  *  matrix  of  values 


INTRODUCTION 


The  effective  use  of  composites  for  high  pressure  cylinders  requires  that, 
the  thick-wall  cylinder  equations  be  rewritten  to  include  cylindrically 
orthotropic  material  properties.  This  process  adds  several  complicating  fac¬ 
tors  to  the  analysis.  First,  is  the  addition  of  two  more  moduli  of  elasticity 
and  two  more  Poisson's  ratios  to  the  material  properties.  Second,  is  that  axial 
stress  in  the  cylinder  is  no  longer  a  constant  in  elastic  analysis.  In  the  new 
case  an  axial  force  will  produce  a  correction  to  all  three  principal  stresses 
which  are  a  function  of  radius. 

Several  factors  which  simplify  the  analysis  of  cylinders  are  common  to  both 
the  isotropic  and  orthotropic  analysis.  The  principal  stress  directions 
correspond  to  the  principal  geometry  directions  of  both  the  cylinder  and  the 
applied  pressures.  This  fact  eliminates  shear  effects  from  the  analysis.  The 
next  simplification  is  that  the  geometry  can  be  nond i mens ionali zed  by  dividing 
by  the  internal  radius  and  it  is  specified  by  the  wall  ratio  (W  *  b/a).  The  two 
loads  are  internal  and  external  pressure  which  are  also  nond i mens i ona 1 .  In 
fact,  stress  can  be  specified  as  a  function  of  either  applied  pressure  when  the 
ratio  of  internal  to  external  pressure  is  a  constant. 

The  form  of  the  equations  used  was  proposed  by  Voyiadjis,  Kiousis,  and 
Hartley  (ref  1)  in  a  1985  paper  dealing  with  the  residual  stresses  in  metal 
cylinders  that  have  been  subjected  to  large  plastic  deformation.  However,  it 
has  been  necessary  to  go  back  to  the  original  work  of  Lekhnitskii  (ref  2)  and 


*G.  Z.  Voyiadjis,  P.  D.  Kiousis,  and  C.  S.  Hartley,  "Analysis  of  Residual 
Stresses  in  Cylindrically  Anisotropic  Materials,"  Experimental  Mechanics,  Vol . 
25,  No.  2,  June  1985,  pp.  145-147. 

*S.  G.  Lekhnitskii,  Theory  of  Elasticity  of  an  Anisotropic  Elastic  Body. 
Holden-Oay  Inc.,  San  Francisco,  1963. 


develop  the  equations  more  fully.  In  the  metallic  problem  the  difference  in  the 
engineering  properties  is  not  very  great  and  some  simplifying  assumptions  are 
possible.  In  contrast,  the  composite  materials  show  a  very  large  difference  in 
properties  with  direction  and  the  full  set  of  equations  is  necessary.  The 
assumptions  that  will  be  used  in  this  report  are  that  of  a  linear  elastic 
cylindrically  orthotropic  material. 

GEOMETRY 

The  problem  is  that  of  a  long  cylinder  (Figure  1)  with  flat  ends.  The 
cylinder  has  a  constant  inner  radius  'a'  and  outer  radius  *b'.  The  ends  are 
fully  constrained  in  the  axial  direction  or  the  classic  plane-strain  condition. 
The  material  is  assumed  to  be  cylindrically  orthotropic  with  its  principal  axes 
coincident  with  the  cylinder.  This  assumption  provides  for  constant  engineering 
properties  in  each  of  the  three  coordinate  directions:  radial,  circumferential, 
and  axial,  while  this  report  will  only  use  the  internal  pressure  load,  the 
equations  are  written  for  constant  pressure  loads  on  either  the  internal  or 
external  surfaces.  Also,  this  report  will  only  use  the  plane-strain  end  con¬ 
dition,  but  will  quote  the  equations  for  axial  loads  necessary  to  calculate 
other  end  load  conditions. 

Thick-Wall  Equations 

For  this  work,  Hooke's  Law  is  used  in  the  following  form: 

[e]  -  [A][a] 

This  must  be  modified  for  use  in  the  cylindrical  problem  and  a  new  material 
compliance  matrix  is  generated  in  using  the  transformation  equation 

Ai3A3i 


In  general,  this  new  matrix  must  be  transforn »d  to  supply  the  proper  values  as 
the  function  of  the  angular  position  in  the  cylinder.  However,  in  the  case  of  a 
cylindrically  orthotropic  material,  this  matrix  has  constant  values. 

In  order  to  solve  thick-wall  cylinder  problems  for  a  variety  of  end  load 
conditions,  i.e.,  the  generalized  plane-strain,  it  is  necessary  to  have  two  sets 
of  equations.  The  first  set  of  four  equations  is  the  plane-strain  solution  for 
a  cylinder  under  a  combination  of  internal  and  external  pressure.  The  second 
set  is  for  the  cylinder  with  an  axial  load  applied.  With  these  equations  the 
stresses  can  be  calculated  for  any  end  load  condition  by  superposition.  Both 
sets  of  equations  contain  a  single  equation  for  stress  in  the  radial,  tangen¬ 
tial,  and  axial  directions.  However,  each  set  contains  a  fourth  equation.  In 
the  pressure  load  set  this  equation  calculates  the  total  axial  load  generated  by 
the  plane-strain  condition.  In  contrast,  in  the  axial  load  set  the  fourth 
equation  calculates  a  rather  complicated  constant  used  in  the  other  three 
equations.  The  eight  equations  were  taken  from  Reference  2  and  corrected  for 
several  minor  typographical  errors,  and  all  terms  were  eliminated  which  are  zero 
for  orthotropic  materials.  The  resulting  equations  are  given  below  at  the  'a' 
equation  in  a  pair  (la,  2a,... 8a).  The  ’b'  equation  is  the  same  equation  in  the 


form  suggested  by  Reference  1.  The  conversion  is  done  in  two  simple  steps. 

b2k  t>2kfc2k 


bk-lfek+l  bk"lbk+lbk*lbk+l 


Next,  make  the  following  substitution: 


*G.  Z.  Voyiadjis,  P.  0.  Kiousi.,,  and  C.  S.  Hartley,  "Analysis  of  Residual 
Stresses  in  Cylindrically  Anisotropic  Materials,"  Experimental  Mechanics,  Voi 
25,  No.  2,  June  1985,  pp.  145-147. 

^S.  G.  Lekhnitskii,  Theory  of  Elasticity  of  an  Anisotropic  Elastic  Body, 
Holden-Oay,  Inc.,  San  Francisco,  1963. 
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The  equations  for  pressure  load  are: 


pak+l  ”  qbk+l  qak“l  “  pbk“l 

Op  »  -- — - — —  pk-1  +  _______ — ____  ak+lbk+lr-k-l 

p2k  _  a2k  b2k  .  a2k 


or  = 


pC0k+1  -  q 

d-c02k) 


r  k-1 

(-)  ♦ 
b 


qC0k_l  -pi  b  k+1 

- !cok+1(-) 

(1-C02k)  J  r 


pak+l  _  qbk+l  qak'l  -  ppk-1 

Oq  * - - — - - krk-1  - - kak+lbk+lr-k-l 

b2k  _  a2k  b2k  _  a2k 


o0  = 


pC0k+1  -  q 
(1-C02k) 

’2  *  *  £ 


r  k-1 

k(-) 

b 


qc0! 


k-1  _ 


(1-C02k) 


.  ,  b  k+1 
kC0k+l(-) 
r 


pak+l  _  qpk+1 

"b2k  _~a2k~" 


(Ai3+kA23)rk_l 


qak-l  -  ppk-1 

♦  .  ak+lbk+1(A13-kA23)r-k-l 


b2k  _  a2k 


'33 


[  (pCQk+l-q)  ] 
(1-C02k) 


(A13+kA23>(-)k‘l  * 
b 


qC0k_l  -  p 
(1-C02k) 


(Ai3-kA23)Cok+1(-) 

r 


k+1 


PP 


2n 

a33 (b2k-a2k) 


(qpk+1 -pak+l ) {bk+l-ak*l ) 


a13  +  kA23 


1  +  k 


+  (qak“l-pbk-^ ) (bk”^-ak“l )a*b*  • 


a13  ~  kA23 
1  -  k  ” 


(la) 


(lb) 


(2a) 


(2b) 


(3a) 


(3b) 


(4a) 


2rt  ...  A13  +  kA23 

PP  -  .  b*(q-PC0k+l)(l-C0k+l)  . 

A33(l-C02k)  1  +  k 


+  a* (qC0k-l-P) (1-C0k“l ) 


a13  -  kA23 
1  -  k 


The  equations  for  axial  load  are: 

Ph  bk+l  -  ak+l  bk-l  -  ak_l 

<7p  S  —  (l - rk_1 - - - ak+lbk+lr-k_l ) 

T  b2k  .  a2k  b2k  _  a2k 


Ph  (l~C0k+l)  r  k-1  (1-C0k-1)  b 

"  1 . (-)  . . c0k+l(-)k+1 

T  (1-C02k)  b  (1-C02k)  r 


Ph  bk+l  -  ak+l  ,  bk-l  -  ak~l 

oq  -  —  (1  - - krk_1  + - kak+lbk+lr"k_l ) 

T  b2k  .  a2k  b2k  _  a2k 


Ph  (1-C0k+1)  r  k-1  (1-C0k“l)  b 

<jQ  *  --  (1 - -  k(-)  ♦ - kC0k+l(-)k+l) 

T  <1-C02k)  b  (1-C02k)  r 


Ph  bk+l  -  ak+1 

'A33  l  A'3  *  -  'lir/li*  <A13*kA23>rk'1 


bk- 1  -  ak”l 
b2k  .  a2k 


ak+lbk+l (Al3_kA23)r-k-l 


Ph  (1-C0k"l)  r 

A13  +  A23 . .  (A13+kA23)  (-) 

TA33  (1-C02k)  b 

{1-Cok'1)  ,  ,  b  k+1 

. I"  (A13-kA23)C0k+1(-)k  1 

1  -  Cn2k  r 


RESULTS 

The  new  solution  has  been  used  to  produce  three  plots.  In  each  of  these 
plots  a  parameter  is  plotted  versus  a  wall  ratio  range  of  1.0  to  2.5  and  is 
calculated  for  the  plane-strain  end  condition.  The  three  parameters  are  solu¬ 
tion  error,  growth  potential,  and  failure  pressure  for  failure  in  simple 
tension.  The  hoop  stress  (circumferential)  at  the  bore  is  used  in  all  cases 
because  this  is  the  critical  stress  in  thick  cylinders.  Solution  error  is  the 


percent  error  when  compared  with  the  standard  strength  of  materials  solution  for 
thin  cylinders.  Growth  potential  is  the  percent  reduction  in  hoop  stress  when 
the  wall  is  increased  in  thickness  by  one  percent.  The  failure  pressure  assumes 
that  the  tube  will  fail  in  hoop  tension  at  a  stress  of  100,000  pressure  units. 
There  are  five  curves  on  each  plot  for  five  values  of  the  orthotropic  material 
parameter  D.  0  is  the  ratio  defined  by  dividing  the  engineering  modulus  in  the 
hoop  direction  by  the  modulus  in  the  radial  direction.  For  0  equal  to  one,  the 
tube  material  is  isotropic  and  for  0  equal  to  128,  the  material  is  highly  ortho¬ 
tropic. 

DISCUSSION 

All  three  of  these  plots  show  the  potential  for  serious  problems  with 
thick-wall  composite  tubes.  The  problems  result  from  the  low  radial  stiffness 
of  the  tubes  relative  to  the  higher  hoop  stiffness.  It  appears  that  the  ratio 
of  hoop  stiffness  divided  by  radial  stiffness  (0  =  Eh/Er)  is  an  important  design 
parameter.  As  this  ratio  becomes  larger,  the  tube  has  diminished  ability  to 
transfer  load  in  the  radial  direction  to  the  outer  fibers.  Therefore,  the  outer 
fibers  may  not  have  the  ability  to  carry  load  as  well  as  the  material  in  the 
same  position  in  an  isotropic  material. 

CONCLUSION 

The  behavior  of  an  orthotropic  high  pressure  cylinder  may  be  very  different 
from  that  of  an  isotropic  cylinder.  Therefore,  any  composite  design  must 
proceed  with  extreme  care.  However,  these  equations  must  first  be  verified 
(Appendix  8). 
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PLANE  STRAIN 


APPENDIX  A 


SOME  RESULTS  ON  ORTHOTROPIC  HIGH  PRESSURE  CYLINDERS 

The  following  four  subprograms  calculate  stresses  in  a  cyl indrical ly 
orthotropic  cylinder: 

ORCYMA  -  is  a  function  that  sets  up  the  material  constants  and  returns  a 
single  parameter  for  use  in  output  identification.  (This  must  be  called  with 
each  change  of  material.) 

ORCYFO  -  is  a  subroutine  that  calculates  the  axial  force  due  to  the  plane- 
strain  end  condition. 

ORCYPR  -  is  a  subroutine  that  calculates  stresses  due  to: 

(a)  internal  pressure 

(b)  external  pressure 

(c)  twisting  about  the  axis 

ORCYTE  -  is  a  subroutine  that  calculates  stresses  due  to  an  axial  force. 


non 


ORCYMA 


C 
C 
C 

C  INPUT 

C  IMAT  =  POINTER  TO  THE  CURRENT  MATERIAL  COMPLIANCE  MATRIX 
C  STRESS  =  PLANE  STRESS  ANALYSIS  ?  (LOGICAL) 

C 

C  IN  THE  /MATERL/  COMMON  BLOCK 

C  AA  =  A  SET  OF  MATERIAL  COMPLIANCE  PROPERTIES 

C 

C  OUTPUT 

C  ORCYMA  =  A  SIMPLIFIED  MATERIAL  CONSTANT  FOR  PRINTER  OUTPUT  ONLY 
C 

C  IN  THE  /MATERL/  COMMON  BLOCK 

C  THE  MATERIAL  CONSTANTS 

C  K  -  K2  -  KP1  -  KN1  -  APKA  -  AMKA  -  A33  -  AI3P23  -  H 
C 

FUNCTION  ORCYMA( IMAT, STRESS) 

IMPLICIT  REAL*8(A-H,K-Z) 

LOGICAL  STRESS 

COMMON  /MATERL/K , K2 , KP1 , KM1 , APKA , AMKA , A33 , A13P23 , H , AA ( 6 , 6 , 1 5 ) 

SET  UP  THE  MATERIAL  CONSTANTS 

All=AA(l, 1, IMAT) 

A22=AA(2,2, IMAT) 

A13=AA(1, 3, IMAT) 

A23= AA( 2,3, IMAT) 

A33=AA( 3,3, IMAT) 

I F ( . NOT. STRESS) THEN 
B11=A11-A13*A13/A33 
B22=A22-A23*A23/A33 
ELSE 
B  1 1  =  A 1 1 
B22=A22 
END  I F 
D=B 1 1/B22 
K=DSQRT( D ) 

KM1=K-1 . 0  DO 
KP 1  =  K+ 1 . 0  DO 
K2-K*2 . 0D0 
APKA=A13+K*A23 
AMKA=A13-K*A23 
A13P23= A13+A23 
H=0 . 0  DO 

IF(BII.NE.B22)H=(A23-AI3)/(BI1-B22) 

ORCYMA=D 
RETURN 
END 
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ORCYFO 


C 
C 

c 

C  INPUT 

C  A  =  INNER  RADIUS 

C  B  =  OUTER  RADIUS 

C  PIN  =  INNER  PRESSURE 

C  POUT  =  OUTER  PRESSURE 

C 

C  OUTPUT 

C  PP  =  PLANE  STRAIN  FORCE 

C 

SUBROUTINE  ORCYFO ( A, B , PIN, POUT, PP) 

IMPLICIT  REAL*8(A-H,K-Z) 

COMMON  /MATERL/K , K2 , KP1 , KM1 , APKA , AMKA, A33 , A13P23 , H 
DATA  PI/3.14159259  DO/ 

CALCULATE  THE  TOTAL  AXIAL  FORCE  FOR  PLANE  STRAIN 
C0=A/B 

PP=B*B*(POUT-PIN*CO**KP1)*( 1 . 0D0-C0»*KP1 ) *APKA/ ( 1 . ODO+K) 
IF(K.NE.l.ODO) THEN 

PP=PP+A*A*(POUT*CO**KM1-PIN)*( 1. 0D0-C0**KMl ) *AMKA/ ( 1. ODO-K) 
ENDIF 

PP=(2.ODO*PI/(A33*(1.0DO-CO**K2)) ) *PP 
RETURN 
END 


> 
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noon 


ORCYPR 


C 
C 
C 

C  INPUT 

C  A  =  INNER  RADIUS 

C  B  =  OUTER  RADIUS 

C  R  =  LOCAL  RADIUS 

C  PIN  =  INNER  PRESSURE 

C  POUT  =  OUTER  PRESSURE 

C  M  =  TWISTING  MOMENT 

C 

C  OUTPUT 

C  SIG  =  PRINCIPAL  STRESSES 

SUBROUTINE  ORC YPR ( A , B , R , PIN , POUT , M , S IG ) 

IMPLICIT  REAL*8( A-H, K-Z) 

COMMON  /MATERL/K , K2 , KP1 , KM1 , APKA , AMKA, A33 , A13P23 , H 
DIMENSION  SIG(4) 

DATA  PI/3.14159259  DO/ 

CALCULATE  THE  STRESSES  FROM  INTERNAL  PRESSURE,  EXTERNAL  PRESSURE 
AND  TORSION  ABOUT  THE  AXIS 

C0=A/B 
C  RADIAL  STRESS 

SIG( 1)  =( (PIN*CO**KPl-POUT)/( 1.0D0-C0**K2) )*( (R/B)**KM1) 

A  +( (POUT*C0**KMl-PIN)/(1.0D0-C0**K2) ) *C0**KP1 * ( B/R ) **KP 1 

C  TANGENTIAL  STRESS 

S IG ( 2 )  =( (PIN*C0**KPl-POUT)/(1.0D0-C0**K2) )*K*( (R/B)**KM1) 

A  -( (POUT*CO**KMl-PIN)/( 1 . 0D0-C0**K2) ) *K*C0**KP 1  * ( B/R ) **K P 1 

C  AXIAL  STRESS 

SIG(3) =( (PIN*C0**KP 1-POUT) /( 1. 0D0-C0**K2) )*APKA*(R/B) **KM1 
SIG(3)=SIG(3)+( (POUT*CO**KM1-PIN)/(1. 0D0-C0**K2) ) *AMKA*C0**KP1 
A  *(B/R)**KP1 
SIG(3)=-SIG(3)/A33 
C  SHEAR  STRESS 

SIG(4)=2.0D0*M*R/(PI*B**4*(1. 0D0-C0**4) ) 

RETURN 

END 
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c 

C  -  ORCYTE  - 
C 

C  INPUT 

C  A  =  INNER  RADIUS 

C  B  =  OUTER  RADIUS 

C  R  =  LOCAL  RADIUS 

C  P  =  TOTAL  APPLIED  FORCE 

C 

C  OUTPUT 

C  SIG  =  PRINCIPAL  STRESSES 
C 

SUBROUTINE  ORCYTE ( A , B , R , P , S IG ) 

IMPLICIT  REAL*8(A-H, K-Z) 

COMMON  /MATERL/K, K2 , KP1 , KM1 , APKA , AMKA , A33 , A13P23 , H 
DIMENSION  S IG ( 3 ) 

DATA  PI/3.14159259  DO/ 

DATA  AOLD,BOLD,KOLD(A3LD/4*O.ODO/ 

C 

C  CALCULATE  THE  STRESSES  FROM  AXIAL  TENSION 
C 

I F ( K . EQ . 1 . 0D0 ) THEN 
S IG ( 1 ) =0 . ODO 
S IG ( 2 ) =0 . 0D0 

SIG(3)=P*1.0D0/(PI*(B*B-A*A) ) 

GO  TO  10 
END  IF 
CO=A/B 

IF ( A . NE . AOLD . OR . B . NE . BOLD . OR . K . NE . HOLD . OR . A33 . NE . A3  LD ) THEN 
T=(B*B-A*A)*A13P23/2.0D0 

T=T-B*B*( ( (1.0D0-C0**KP1 } **2/(1. 0D0-C0**K2) )*( APKA/ (K+l. ODO) )  ) 

T  =  T-A*A*( ( ( 1.0D0-C0**KM1)**2/(1.0D0-C0**K2) )*( AMKA/ ( K- 1 . ODO )  )  ) 

T=PI*(B*B-A*A)-(2. 0D0*PI*H/A33)*T 

AOLD= A 

BOLD=B 

KOLD-K 

A3LD=A33 

ENDIF 


S IG ( 1 ) = 1 . ODO- ( ( 1.0D0-C0**KP1)/(1.0D0-C0**K2) )*(R/B)**KM1 
SIG(l)=SIG(l)-( (1. 0D0-C0**KM1 )/(1.0DO-C0**K2) )*C0**KP1*(B/R)**KP1 
SIG(1)=(H/T)*SIG(1)*P 

SIG(2)=1.0D0-((l.0D0-C0**KPl)/(1.0D0-C0**K2))*Kt(R/B)**KMl 
SIG(2)=SIG(2)+( ( 1.0D0-Cn**KMl)/( 1 . 0D0-C0**K2 ) )*K*C0**KP1* 

A  ( B/R ) **KP 1 
SIG(2)=(H/T)*SIG(2)*P 

SIG(3)=Al3P23-( ( 1.0D0-C0#*KP1)/(1.0D0-C0**K2) ) *APKA* ( R/B ) **KM1 
SIG(3)=SIG(3)-( (1.0D0-C0**KM1)/(1.0D0-C0**K2) ) *AMKA«C0**KP1 * 

A  ( B/R) **KP1 

SIG(3)=(1.0D0/T-(H/(T*A33) )*SIG(3) )*P 
10  RETURN 
END 
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APPENDIX  B 


SOME  RESULTS  ON  ORTHOTROPIC  HIGH  PRESSURE  CYLINDERS 
SOLUTION  VERIFICATION 

This  solution  was  verified  by  comparison  with  a  ten-element  finite  element 
solution  using  the  ABAQUS  (4.5-175)  code.  ABAQUS  is  a  general  nonlinear  program 
supplied  by: 

Hibbitt,  Karlsson,  and  Sorensen,  Inc. 

101  Medway  Street 
Providence,  RI  02906 
(401)  861-0820 

For  this  problem  only  the  linear  portion  of  ABAQUS  was  used  in  combination  with 
an  eight  node  axisymmetric  element  (CAX8)  and  a  linear  elastic  orthotropic 
material  definition.  It  was  necessary  to  use  only  a  single  row  of  ten  elements 
for  all  solutions  and  the  stresses  were  recovered  at  the  nodal  points.  The 
proper  use  of  nodal  point  constraints  and  constraint  equations  insured  that  the 
radial  faces  of  the  model  remained  flat  and  parallel. 

The  material  properties  were  that  of  a  graphite  epoxy  wrapped  composite 
cylinder.  The  material  was  used  in  a  major  program  at  NASA  (ref  B-l).  For  con¬ 
venience,  the  theoretical  analysis  program  was  written  to  use  the  same  material 
definition  record  as  ABAQUS.  This  also  insured  that  both  programs  used  exactly 
the  same  material  information. 

A  total  of  18  ABAQUS  solutions  were  generated  for  a  selection  of  three  dif¬ 
ferent  wall  ratios,  four  different  loading  conditions  and  two  fiber  wrap  angles 
(material  properties).  Four  solutions  have  been  selected  and  the  results  are 
shown  in  Tables  B-l  through  8-4.  The  four  cases  were  selected  for  no  particular 
reason  except  that  each  of  the  four  loads  is  represented. 

8_1V.  Verderaime,  "Development  of  In  Situ  Stiffness  Properties  For  Shuttle 
Booster  Filament  Wound  Case,"  NASA  Technical  Paper  2377,  George  C.  Marshall 
Space  Flight  Center,  Huntsville,  AL,  August  1984. 


TABLE  «-l 


0B  THDTROP  IC  HIGH  PRESSURE  CYLINDER  THEORY 

SOLUTION  COMPARISON  WITH  ABACUS  USING  10  CAX8  ELEMENTS 

FOR  4  GRAPHITE  EPOXY  OPEN  END  CYLINDER 

WALL  RATIO  =  1.50  -  wRAP  ANGLE  =  0.0  DEGREES 

APPLIEO  L040S 
INTERNAL  PRESSURE  *  1.0 
EXTERNAL  PRESSURE  =  0.0 
AVERAGE  AXIAL  STRESS  *  0.0 


MATERIAL  STIFFNESS  MATRIX 


XX 

ZZ 

YY 

XZ 

XY 

ZY 

XX 

3.293E+9 

3. 586F+9 

3.791E+9  - 

O.OOOE+O 

0  . 0  00  £  +  J 

0 . OOOE  +  0 

LL 

3.536E+9 

3.298F+9 

3.79LF  +  9  - 

o.oooe+o 

O.OOOE+O 

O.OOOE+O 

YY 

3. 79 1 E  +  9 

3.791E+9 

122.4E+9  - 

O.OGQE+O 

o.oooe+o 

0 . 0  0  JE  +0 

XZ 

O.OOOE+O 

0  •  0  00  £  ♦  0 

O.JOOF+0  - 

2 • 550E+9 

O.OOOE+O 

O.OOOE+O 

XY 

O.OOOE+O 

0.0006+* 

O.OOOE+O  - 

o.oooe+o 

6.54QE+9 

O.OOOE+O 

LX 

O.OOOE+O 

O.OCOS+O 

0*  OOOE +  0  - 

o.oooe+o 

u.OOOE+O 

6 . 540F  +c 

CALCULATED  STRESSES 


P  ADI  AL 

STRESS 

HOD  P 

STRESS 

AX  IAL 

STRESS 

RADIUS 

THEORY 

ABACUS 

THEORY 

ABACUS 

THEORY 

A  -nous 

I.  000 

-l.  000 

-.994 

4.330 

4.33 

-.253 

-  c  25  0 

1.050 

-.  76  3 

-.  762 

3.474 

3.48 

-.  170 

157 

l.  100 

-.  590 

-.586 

2.33  3 

2.83 

-.  106 

-.104 

1.150 

-.452 

-.449 

2.347 

2.35 

-o  056 

-c  05  5 

1.200 

-o  344 

-.342 

1.978 

1.98 

-  c  0 1 6 

—  .018 

l.  250 

—  c  25  7 

-.255 

1.695 

1. 70 

0  o  0  l  5 

)  o  0 1  c 

I.  300 

-.136 

-.  135 

1.478 

l  •  48 

0.041 

0.042 

1.350 

-.  128 

-.127 

1.313 

1.31 

0.  053 

0 . 5  '•<3  J 

1.400 

-.078 

-.  078 

l.  13  3 

1.1  $ 

Oo  08  2 

j  .  0  J  2 

1.450 

-.036 

-.  036 

1.094 

1.09 

0.0*3 

0,C9a 

1.500 

0.000 

0.000 

l.  02  6 

1.03 

Oo  112 

0.112 
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TABLE  3-2 

ORTHOTRGPIC  HIGH  PRESSURE  CYLINDER  THEORY 

SOLUTION  COMPARISON  WITH  ABACUS  USING  10  CAX3  ELEMENTS 

FOR  A  GRAPHITE  EPOXY  CYLINDER  WITH  END  LOADS 


WALL  0  A  T I Q  =  1.75  -  , 

A  P°L I FO  LOADS 
INTERNAL  PRESSURE  =  0.0 
EXTERNAL  PRESSURE  *  0.0 
AVERAGE  AXIAL  STRESS  =  1.0 


MATERIAL  STIFFNESS  MATRIX 

XX  2  Z  YY 

XX  3.293E+9  3.610E+9  3.767E+9 

ZZ  3.610E+9  11.63E  +  P  13.31E+9 

YY  3.767E+9  13.31E  +  9  99.05E+9 


WRAP  ANGLES  =  20. 0»  — 2  0. 0  DEGREES 


xz 

O.OOOE+3 

O.OOOE+O 

O.OOOE+O 


XY 

O.QOOE+O 
O.OOOE+O 
0  .000E+ 0 


ZY 

0. DOOR +0 
■JeGw' JEO 

:.o:  of  +  o 


xz  O.OOOE+O 
XY  O.OOOE+O 
ZY  O.QOOE+O 


0 .000  E  +0 

o.oooe+o 

O.OOOE+O 


O.OOOE+O 
O.OOOE+O 
0. 30Q£  +  3 


CALCULATED  STRESSES 


8 AOI US 
1.000 
1.075 
1.  150 
1.225 
1.300 
l.  375 
1.450 
1.525 
l.fcOO 
1.675 
1.750 


RADIAL 
THEORY 
0.000 
0.043 
0.  06  7 
0.  078 
0.031 
0.077 
0.06  7 
0.056 
0.040 
0.021 
0.000 


STRESS 
ABAOUS 
0.002 
J.044 
0.068 
0.  079 
0.082 
0.078 
0.069 
0.056 
0.040 
0.021 
0.000 


HOOP 
THEORY 
Oc  733 
0.  503 
0.328 
0.137 
0.067 
-.  040 
140 
-.  236 
-.  332 
-.42  8 
-.526 


3.Q17E+9 

O.OOOE+O 

O.OOOE+O 


STRESS 
ABAOUS 
0.  734 
0.  504 
0.329 
0.188 
0.068 
0.040 
0.140 
0.236 
0.332 
0.  42  3 
3.526 


O.COOE+O 
6 .373E+- 
O.OOOE+O 


0  .OPOE+O 
0.1  OCF  +  3 

l6.8ac  +  '7 


AXIAL  STRESS 
THEORY  ABAOjS 


1.C73 
1.06  3 
1.053 
1.040 
l.  02  6 
1.0L1 
0.995 
0.973 
U .  ~  6  0 
0.941 

0.521 


1 . ;  8 

1.3  7 
1.05 
1.04 
1.  33 
l  c  31 

0.7  7  C 
0  .  c  7  w 

0.*  6C 
3.9  +  1 
0.92  1 
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TABLE  B-3 


ORTHCTROPIC  HIGH  PRESSURE  CYLINDER  THEORY 


SOLUTION  COMPARISON  WITH  ABACUS  USING  10  CAX8  ELEMENTS 
FOR  A  GRAPHITE  EPOXY  PLANE  STRAIN  CYLINDER 


WALL  RATIO  *  2.00  -  WRAP  ANGLE—  0.0  DEGREES 

APPLIED  LOAOS 
INTERNAL  PRESSURE  =  0.0 
EXTERNAL  PRESSURE  =  1.0 
AVERAGE  AXIAL  STRESS  *  0.0 


MATERIAL  STIFFNESS  MATRIX 

XX  ZZ  YY 

XX  3.298E+9  3.586E+9  3.79l=+9 

ZZ  3. 536E +9  3.298E+9  3.7912+9 

YY  3.791E  +  9  3.791E+9  122.4E+9 


XL 

O.OOOE+O 

O.OOOE+O 

O.OOOE+O 


XY 

O.OOOE+O 
O.OJOE+  ’) 
O.OOOE+O 


ZY 

0. OOOt+O 
J . jOOE +  w 
n„ 100E+0 


XL 

O.OOOE+O 

O.OOOE+O 

0.3G0E  +  0  - 

2.550E+9 

O.OOJE+O 

0.00 jE+O 

XY 

O.OOOE+O 

0.0006+0 

0.  OOOE  +  O  - 

O.OOOE+O 

6.5406+3 

0.0002+0 

zx 

O.OOOE+O 

O.OOOE+O 

O.OOOE+O  - 

O.OOoE+O 

0.0006+0 

6.6^0E+9 

calculated 

i  STRESSES 

RAOIAL  STRESS 

HQUP 

STRESS 

AXIAL 

STRESS 

RADIUS 

THEORY 

ABACUS 

THEORY 

ABACUS 

THEORY 

ABACUS 

1.000 

0.000 

-.002 

-1.078 

-1.08 

-.019 

— .  u  ?  o 

1.100 

-.095 

-.097 

-1.046 

-1.05 

-.059 

-.D  60 

1.200 

-.177 

-.173 

-1.127 

-1.13 

-.09  5 

-.09  6 

1.300 

-.256 

-.256 

-1.287 

-1.29 

-.132 

-.1  37 

1.400 

-.337 

-.337 

-1.  507 

-l.  51 

-.170 

-.170 

l.  500 

-.424 

-.424 

-1.760 

-1.73 

-.211 

-.211 

1.600 

-.519 

-.513 

-2.103 

-2.  10 

-.257 

-.257 

1.700 

-.623 

-.622 

-2.474 

-2.47 

-.  308 

-.3  OR 

1.300 

-.737 

-.  736 

-2.892 

-2.39 

-.364 

- .  3  64 

1.900 

-.362 

-.861 

-3. 360 

-3.36 

-  •  *?o 

-.42  * 

2.  000 

1.000 

0.999 

-3.878 

-3.38 

-.49  3 

-.45  3 
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orthutrqpic  high  pressure  cylinoep  theory 

SOLUTION  COMPARISON  WITH  ABAQUS  USING  10  CAX8  ELEMENTS 
C0P  A  GRAPHITE  EPOXY  PLANE  STRAIN  CYLINDER 


WALL  RATIO  ■  2.00  -  WRAP  ANGLE  =  0.0  DEGREES 

ALLIED  LOADS 
INTERNAL  PRESSURE  =  1.0 
EXTERNAL  PRESSURE  *  1.0 
AVERAGE  AXIAL  STRESS  =  0,0 


MATERIAL  STIFFNESS  MATRIX 


XX 

XX  8.298E+9 
LI  3.5S6E+9 
YY  3 • 7° 1 E+9 


ll 

3.536E+9 
3 .298E+9 
3 . 791 E ♦9 


YY 

3. 791 E  +9 
3.791E+9 
122.4E+9 


XZ 

O.OOOE+O 

O.OOOE+O 

O.OOOE+O 


XY 

O.OOOE+O 
0 . OOOE  +  J 
O.OOOE+O 


ZY 

).  OOOE+J 
G  .  0  0  0  E  +  J 
0.  900F+0 


XZ  O.OOOE+O  O.OOOE+O 
XY  J.OOOE+O  O.OOOE+O 
ZX  O.OOOE+O  O.OOOE+O 


0 . D  0  0  E +0 
O.OOOE+O 
O.OOOE+O 


2. 550E+9 
O.OOOE+O 
O.OOOE+O 


o.Dooe+o 

6.540E+9 

O.OOOE+O 


0  •  0  '■*  0  E  +  J 
0 .  U  C  0  E  +  0 

6  .  54-JF+c 


CALCULATED  STRESSES 

RADIAL  STRESS 


R  AO  I  US 

THEORY 

ABAQUS 

l.  JO'J 

1 . 00  0 

-.981 

l.  LOO 

-.723 

-.709 

I.  200 

-.  585 

-.577 

1.  300 

528 

-.523 

1.400 

-.522 

-.519 

1.500 

-.550 

-.548 

l.  600 

-.603 

-.602 

I.  700 

-.678 

-.676 

1.300 

-.  769 

-.768 

1.900 

-.  877 

-.876 

2.000 

1.000 

0.999 

HOOP 

STRESS 

theory 

ABAQ'JS 

2.801 

2.  81 

1.412 

1.42 

0.501 

0.505 

-.  163 

-.  lol 

-.701 

-.699 

-l.  179 

-1.18 

-1.635 

-1.63 

-2.093 

-2.09 

-2. 563 

-2.57 

-3.071 

-3.07 

-3.609 

-3.61 

AXIAL 

STRESS 

theory 

1  L>A  5u  S 

-.374 

—  •  3  c  £ 

-.23  I 

—  .27'; 

-.209 

-.2?> 

-.227 

-.225 

-.234 

-.2  32 

-.254 

-c  25  ^ 

-.23  5 

-.’H 

-.325 

-.324 

-.372 

-.371 

-.427 

—  o  ff  2  ^ 

-.  4-»>? 

-.43  V 
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TECHNICAL  REPORT  INTERNAL  DISTRIBUTION  LIST 

NO,  OF 
COPIES 

CHIEF,  DEVELOPMENT  ENGINEERING  BRANCH 

ATTN:  SMCAR-CCB-D  1 

-OA  i 

-OC  3 

-OM  i 

-DP  i 

-DR  i 

-DS  (SYSTEMS)  1 

CHIEF,  ENGINEERING  SUPPORT  BRANCH 

ATTN:  SMCAR-CCB-S  ! 

-SE  1 

CHIEF,  RESEARCH  BRANCH 

ATTN:  SMCAR-CCB-R  2 

-R  (ELLEN  FOGARTY)  ! 

-RA  i 

-RM  i 

-RP  ! 

-RT  1 

TECHNICAL  LIBRARY  5 

ATTN:  SMCAR-CCB-TL 

TECHNICAL  PUBLICATIONS  &  EDITING  UNIT  2 

ATTN:  SMCAR-CCB-TL 

DIRECTOR,  OPERATIONS  DIRECTORATE  1 

ATTN:  SMCWV-OD 

DIRECTOR,  PROCUREMENT  DIRECTORATE  1 

ATTN :  SMCWV-PP 

DIRECTOR,  PRODUCT  ASSURANCE  DIRECTORATE  1 

ATTN :  SMCWV-QA 


NOTE:  PLEASE  NOTIFY  DIRECTOR,  BENET  WEAPONS  LABORATORY,  ATTN:  SMCAR-CCB-TL 
OF  ANY  ADORESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 

ASST  SEC  OF  THE  ARMY 

RESEARCH  AND  DEVELOPMENT 

ATTN:  DEPT  FOR  SCI  ANO  TECH  1 

THE  PENTAGON 

WASHINGTON,  D.C.  20310-0103 
ADMINISTRATOR 

DEFENSE  TECHNICAL  INFO  CENTER 
ATTN:  DTIC-FDAC  12 

CAMERON  STATION 
ALEXANDRIA,  VA  22304-6145 


COMMANDER 
US  ARMY  AROEC 

ATTN:  SMCAR-AEE  1 

SMC AR- AES,  BLOG.  321  1 

SMCAR-AET-O,  BLDG.  351N  1 

SMCAR-CC  1 

SMCAR-CCP-A  1 

SMCAR-FSA  1 

SMCAR-FSM-E  1 

SMCAR-FSS-O,  BLOG.  94  1 

SMCAR-MSI  (STINFO)  2 

PICATINNY  ARSENAL,  NJ  07806-5000 

DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN:  SLCBR-DD-T,  BLDG.  305  1 

ABERDEEN  PROVING  GROUND,  MO  21005-5066 

DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS  ACTV 
ATTN:  AMXSY-MP  1 

ABERDEEN  PROVING  GROUNO,  MO  21005-5071 

COMMANDER 
HQ,  AMCCOM 

ATTN:  AMSMC-IMP-L  1 

ROCK  ISLANO,  IL  61299-6000 


NO.  OF 
COPIES 

COMMANOER 

ROCK  ISLANO  ARSENAL 

ATTN:  SMCRI-ENM  1 

ROCK  ISLAND,  IL  61299-5000 

OIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENGR  ACTV 
ATTN:  AMXIB-P  1 

ROCK  ISLAND,  IL  61299-7260 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMMAND 
ATTN:  AMSTA-DDL  (TECH  LIB)  1 

WARREN,  MI  48397-5000 

COMMANDER 

US  MILITARY  ACADEMY  1 

ATTN:  DEPARTMENT  OF  MECHANICS 
WEST  POINT,  NY  10996-1792 

US  ARMY  MISSILE  COMMANO 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 
REDSTONE  ARSENAL,  AL  35898-5241 

COMMANDER 

US  ARMY  FGN  SCIENCE  ANO  TECH  CTR 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  LABCOM 

MATERIALS  TECHNOLOGY  LAB 

ATTN:  SLCMT-IML  (TECH  LIB)  2 

WATERTOWN,  MA  02172-0001 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT'O) 


NO.  OF 
COPIES 

COMMANOER 

US  ARMY  LABCOM,  ISA 

ATTN:  SLCIS-IM-TL  1 

2800  POWOER  MILL  ROAO 
AOELPHI,  MO  20783-1145 

COMMANOER 

US  ARMY  RESEARCH  OFFICE 

ATTN:  CHIEF,  IPO  1 

P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709-2211 
OIRECTOR 

US  NAVAL  RESEARCH  LAB 

ATTN:  DIR,  MECH  DIV  1 

CODE  26-27  (DOC  LIB)  1 

WASHINGTON,  O.C.  20375 


NO  OF 
COPIES 

COMMANOER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MN  1 

EGLIN  AFB,  FL  32543-5434 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MNG 

EGLIN  AFB,  FL  32542-5000  1 

METALS  AND  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  DIVISION 
505  KING  AVENUE 

COLUMBUS,  OH  43201-2693  1 


NOTE:  PLEASE  NOTIFY  COMMANOER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  WEAPONS  LABORATORY,  SMCAR-CCB-H , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADORESS  CHANGES. 


